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Abstract 

Light-front quantization has important advantages for describing relativistic statisti- 
cal systems, particularly systems for which boost invariance is essential, such as the 
fireball created in a heavy ion collisions. In this paper we develop light-front field 
theory at finite temperature and density with special attention to quantum chromo- 
dynamics. We construct the most general form of the statistical operator allowed by 
the Poincare algebra and show that there are no zero-mode related problems when 
describing phase transitions. We then demonstrate a direct connection between den- 
sities in light-front thermal field theory and the parton distributions measured in 
hard scattering experiments. Our approach thus generalizes the concept of a parton 
distribution to finite temperature. In light-front quantization, the gauge-invariant 
Green's functions of a quark in a medium can be defined in terms of just 2-component 
spinors and have a much simpler spinor structure than the equal-time fermion prop- 
agator. From the Green's function, we introduce the new concept of a light-front 
density matrix, whose matrix elements are related to forward and to off-diagonal par- 
ton distributions. Furthermore, we explain how thermodynamic quantities can be 
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calculated in discretized light-cone quantization, which is applicable at high chemical 
potential and is not plagued by the fermion-doubling problem. 
PACS: ll.lO.Wx, 12.38.Lg, 12.38.Mh, 24.85.+p 

Keywords: Light-Front Quantization, Thermal Field Theory, Generalized Parton 
Distributions 



I. INTRODUCTION 

Dirac's front form of relativistic dynamics Q] has remarkable advantages for relativistic 
problems in high energy and nuclear physics. Most appealing is the simplicity of the vacuum 
(the ground state of the free theory is also the ground state of the full theory) and the 
existence of boost-invariant light cone wavefunctions (see Ref. j2| for a review.) This makes 
light-front quantization a natural candidate for the description of systems for which boost 
invariance is an issue, such as the fireball created in a heavy ion collision or the small- a; 
features of a nuclear wavefunction. It is clearly important to exploit the advantages of light- 



front quantization for thermal field theorj^ 
been done by several authors 




uable work in this direction has already 



In this paper we shall apply light- 



front quantization to statistical physics and investigate the prospects and challenges of this 
approach for quantum chromodynamic systems. 

In the front form, see Fig. Q initial conditions are defined on a light-like hypersurface 
with 7-+ = -)- = Q_ ^\ summary of our notation is given in Appendix El) This sets the 
boundary condition as the light front moves forward. When the theory is quantized, (anti- ) 
commutation relations are defined at fixed light-front time r+, instead of fixed equal-time r°. 
As a consequence, the number of kinematic Poincare generators in the front-form is larger 
than in any other form |2|, namely 7 out of 10 Poincare generators do not depend on the 
interaction. In particular, the generator of Lorentz boosts in the direction is a kinematic 
operator on the light front, but not in the instant form. Therefore, the statistical operator 
in light-front quantization will transform trivially under longitudinal boosts. 

Another advantage of the light-front formulation of thermal field theory is that the charge 
is defined as the integral over the +-component of the current density, Q = J d^rj^{r). This 
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FIG. 1: The three forms of Hamiltonian dynamics jj, [J/. In the front form (middle), initial 
conditions are defined on a light-like surface, and the system is the propagated in light-cone time 
x'^ by the light-cone Hamiltonian P~ . 



allows one to establish a clear connection between densities in thermal field theory and 
the boost-invariant parton distributions measured in hard scattering experiments [1^. No 
such relation exists in the instant form, where charges are integrals over 0-components of 
current densities. Also, parton distributions defined in equal-time quantization are not boost 
invariant, since the generator of Lorentz transformations in direction contains interaction 
terms. As a consequence, in this formalism one might regard low-x gluon saturation as 
a phase transition 



1^ or a critical phenomenon in the sense of statistical mechanics. 
One can even generalize the concept of a parton distribution to finite temperature, which 
is relevant for jet quenching in heavy ion collisions as well as for parton recombination 



We also point out that fermions behave completely differently on the light-front, where 
they are described terms of just 2-component spinors. It turns out that the fermion propa- 
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gator has only one derivative in the numerator, leading to only one pair of doublers on the 
lattice. In addition, the doubler problem is entirely avoided in Discretized Light-Cone Quan- 
tization (DLCQ) [l8|. This technique is also applicable at finite density and will hopefully 
allow one to study baryon-rich QCD matter in large scale numerical calculations. Despite 
recent progress [l9|, conventional lattice QCD is still notoriously complicated at finite baryon 
chemical potential. 

This paper is organized as follows. In the next section, we construct the most general 
form of the statistical operator allowed by the Poincare algebra. Our results are in agreement 
.tH ea.„. and... ffl Q. We a..ue .Ka. ..e .n,p„c.. of .Ke _ doe. „ot lead 
to difficulties when describing phase transitions. In Section IIIIl we define gauge-invariant 
Green's functions for fermions at finite temperature and relate them to generalized parton 
distributions in Section HVl In the summary, we give an outlook on further applications and 
challenges for this approach. 

II. THE STATISTICAL OPERATOR AT FINITE TEMPERATURE AND 

DENSITY 

The form of the statistical operator w at finite temperature and density can be obtained 
from very general considerations. Our result below for w, Eq. ®, is compatible with the 



find... of Ba, is a,wa,s .he exponential of the equal t,.e ene,, P« . the loeal .e.t 

frame of the system. Here, we shall give a more physical derivation of that result, following 
Ref. Q- 

In the instant form, the equal time Liouville theorem [iB^wet = [P^,'Wet]) requires 
that in equilibrium, wet is a function of only those Poincare generators which commute 
with the equal time Hamiltonian P°. In addition, since systems far apart from each other 
must be uncorrelated, the density operator of the combined system has to factorize into the 
density operators of the subsystems. This requirement is similar to the cluster decomposition 
principle that has been studied in the context of the deuteron wavefunction on the light- 
cone in Ref. [2^. Consequently, in equilibrium Itl^wet) must be a linear combination of 
the additive constants of motion, namely the four components of the momentum P'^ and 
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the angular momentum three-vector J. Hence, \n.{wET) = a — PiuyP" — uj • J — Ylil^iQi)-! 
where a is a normahzation constant. The notation is chosen in anticipation of the physical 
meaning of the various coefficients: /3 is the inverse temperature and is the four velocity 
of the system, cf. Ref. In addition, u is the angular velocity at which the body rotates. 
Additional conserved charges Qi are included along with their chemical potentials /i;. We 
stress that the "charges", i.e. the mean values of the operators P^,, J, Qi, are given and that 
the "chemical potentials" {u,y,uj,fii) need to be determined from the conservation laws. In 
quantum mechanics, of course, one can simultaneously specify only charges which commute 
with each other, i.e. only one component of J (namely for rotations around the axis with the 
largest moment of inertia.) For the same reason, one cannot specify all four P^, for systems 
with nonzero angular momentum. 

What changes, when we switch to the light front? First of all, the equal time Liouville 
theorem has to be replaced by its light front analog. Writing the light front density operator 
as 

the (light cone) time {x~^) evolution of the multi particle state {(pi) is generated by the second 
quantized Hamiltonian P~, 

td- 10,) = P- 10,) . (2) 
From this, the light front version of the Liouville theorem immediately follows: 



zd w 



(3) 



Equilibrium (on the light front) is reached, when the commutator on the right hand side 
(rhs) of Eq. ^ vanishes. The six Poincare generators which commute with P~ are the other 
three momentum components and three of the six independent elements of the antisymmetric 
boost-angular momentum tensor 

/ n + Jo kn - J, ~2k. \ 



i^i + Js K2 - Ji -2i^3 
-K1-J2 J3 -K1 + J2 

-K2 + J1 -J3 -K2-J1 

2K3 k, -% k^ + Ji 



(4) 
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(/i, V G {+, 1, 2, — }). Namely the rotations (M^'^ = J3) around the longitudinal direction and 
the two dynamic operators M^' = (M^", M^-), where M^^ = -Ki + J2, M^^ = -K2 - Ji- 
The Kj are boost generators in j direction. Since two of the operators M-^~ have zero 
eigenvalues when acting on all physical states, the light front density operator has the general 
form 

ln(e) = a-P (^a^P'' -uJ.-J^ f^iQi^ ■ (5) 

Thus we have obtained the same Jiittner type distribution as in the instant form. The 
identification of with the four-velocity of the system relative to the observer can be 
derived by calculating the commutators of ln('u;) with the boost generators by means of the 
Poincare algebra. Since u^u'^ = 1, only four of the five parameters (3 = 1/T and u'^ are 
independent. 

We remark that T is the same as the instant-form temperature, but the chemical po- 
tential has a different meaning. On the light-front, densities are given by +-components of 
currents, and not by 0-components. This is essential for the proper generalization of parton 
distributions to finite temperature, since the latter are also defined as +-components jl^. 

Note that one can include rotations in the light front formalism, but one has to make sure 
that the longitudinal axis has the largest inertial moment. Of course, one cannot specify 
and J3 simultaneously. We shall not study systems with finite angular momentum in what 
follows. 

We choose a = as normalization so that the partition function is given hy Z = Tiw. 
Since Z is a. Lorentz scalar, all thermodynamic potentials and the entropy transform as 
scalars, e.g. the Lorentz invariant generalization of the grand-canonical potential (or of the 
free energy in the case of /i = 0) is 

Vt = -T\nZ{y,T,n), (6) 

and the entropy is defined by 

The role of the total energy of the system is now played by the expectation value of UyP^ ^ 

7/ ^ 

U={u,P'') = ^i:i{wP'). (8) 
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As usual, f2 = U — TS — nQ. All known relations between thermodynamic potentials remain 
valid. 

Even though one can formulate light-front statistical mechanics in an arbitrary frame, 
it is convenient to choose the reference frame in which the system as a whole is at rest, so 
that ln(u}) = —(3P^. Since the three forms of Hamiltonian Dynamics are not related by 
Lorentz transformations, there is no frame in which = 0. Therefore, w 7^ exp{—j3LcP ) 
(cf. Ref. 0,3). 

The quantities (5, u^, 00 and fii, have the meaning of Lagrange multipliers that hold the 
mean values of the constants of motion fixed, while entropy is maximized. In an ideal gas 
for example, the maximum entropy is attained for occupation numbers given by Fermi-Dirac 
and Bose-Einstein statistics 



;ntropy 

flfl 



assuming that particles carry charge +1 and antiparticles charge —1. (Note that u^p'^ > 0.) 
Here, p'^ is the 4-momentum of a single particle and the degeneracy factor for different spin 
states is denoted by g. The Lagrange multipliers define the equilibrium conditions for two 
systems. In complete equilibrium with each other, both systems must have the same values 
of temperature, u^, u and /i;, i.e. no internal motion of macroscopic parts of the system is 
possible in equilibrium (at least in the absence of vortex lines In particular, the system 

cannot perform a relative motion with respect to an external heat bath. It is important in 
this respect that the values of these Lagrange multipliers do not depend on the total size of 
the system. Thus, even though one can rewrite the statistical operator in terms of light-front 
momentum fractions xi and the total invariant mass squared (M^) as (let P± = 0± and all 

Qi = o) 

w = exp(-7^ - iM^M"^) (10) 

I 

the quantities 'j^ = u~P'^/{2T) and 7j\/2 = m+/(2TP+) depend on the overall size of the 
system and are therefore unsuitable to define equilibrium conditions. For that reason, some 
care is required when talking about phase dia gram s in the x-M^ plane, which is often done 
in the context of the Color Glass Condensate |lJ]. 
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As usual, the (grand) canonical ensemble allows for fluctuations of the total momentum 
and of the total charge. Since we have constructed w such that it depends only on extensive 
quantities, the relative magnitude of these fluctuations becomes irrelevant for large systems 
such as a neutron star, i.e. 



for each u. Here N is the number of quanta in the system. For a large enough system 
one can therefore use the canonical or the grand canonical ensembles, which are much more 
convenient for calculations than the microcanonical ensemble, even if the values of , J3 
and Qi are strictly conserved. For a large nucleus, > 10^. This means a statistical 
approach is still justified. For a single hadron however, the microcanonical ensemble with 
strict conservation of the constants of motion is more appropriate [2^ . 

The simplicity of the light-front vacuum, usually considered an advantage, seems to bear 
problems as far as phase transitions are concerned. Clearly, thermal field theory on the light- 
front would be rendered useless if it turned out that the zero-mode problem [2^ [2^ would 
have to be solved before any statements about phase transitions could be made. However, the 
fact that the statistical weight of a configuration is maximized for minimal equal-time energy 
rather than for minimal light-front energy has important consequences for the formation of 
condensates at low temperatures. The ground state, i.e. the state the system is in at T = 0, 
is in general different from the light-front vacuum. (Indeed, for massive particles P~ = 
corresponds to infinite P^-) 

Alves, Das and Perez |2| obtained the self-energy to one-loop order in 0^ theory [C = 
dijcf)* d'^ (j) — m^(j)*(f) — A(0*0)^). Replacing —m? with > 0, one immediately obtains spon- 
taneous symmetry breaking from their results, including a massless Goldstone boson. No 
problem arises from l//c~'^-poles. As usual, one of the fields in this model develops a non- 
vanishing expectation value and is shifted correspondingly. One of the redefined fields, the 
(J, acquires a finite mass. The other field is massless and will be referred to as vr here. In 
the high temperature limit, the masses obtained from the one- loop calculation of Ref. [3] are 
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given by the simple expressions 

2 .'2c2(l-TVT,2), T<T. 

= < (12) 
|A(T2-T2), T>T, 

ml = { - , (13) 

|A(T2-T2), T>r, 

with Tc = 3c^/A. These are the standard results 2^. The broken phase has one massive 
sigma and one massless pion. There is a second order phase-transition at T = T^, where all 
masses vanish. In the symmetric phase, all fields have the same mass. 

In addition, the authors of Ref. jl^ reproduced the chiral phase transition in the Nambu- 
Jona-Lasinio Model on the light-front. 

The possibility of spontaneous symmetry breaking without zero-modes has previously 



been suggested in Ref. [26[ for 0^(1+1), but these authors did not study the finite tem- 
perature case. Nevertheless, the picture arising in our work is similar to that of Ref. j^l, 
where P~ is minimized, keeping the total fixed by introducing a Lagrange multiplier. In 
our approach, on the other hand, the equilibrium configuration is found by maximizing the 
entropy for fixed mean values of all four momentum components. 

We conclude that this approach is poised for the study of phase transitions in more 
complicated field theories, such as QCD. 



III. FERMION GREEN'S FUNCTION IN A MEDIUM 

Until now one could get the impression that thermodynamics and statistical physics on 
the light front is identical to the usual instant form approach, except for a trivial change of 
variables. That this is not the case becomes most clear, when one studies fermions on the 
light-front. 

In light-front field theory, the Dirac equations can be written as a set of two coupled 
equations for 2-component spinors, see appendix El Only one of these equations contains a 
time derivative, the other one is a constraint. As a consequence, the entire theory can be 
formulated in terms of 2-component spinors, very much like a non-relativistic theory. It is 
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however important to note, that the 2-component theory does not follow from a local, scalar 
Lagrangian. The theory is non-local, because the action can propagate "instantaneously" 
along the light-line, where the quantization surface touches the light-cone (see Fig. P). This 
often leads to difficulties known as the zero-mode problem. We believe that these difficulties 
can be overcome by introducing time-ordered, retarded and advanced Green's functions 
directly in the Hamiltonian formulation. This is done in Ref. [2^ for the non-relativistic 
case, in which non-local operators are common, since particles interact through potentials. 
Nevertheless, in a diagrammatic expansion of these Green's functions, all 4 momentum 
components are conserved at each vertex. The Green's function is the fundamental object 
of this approach. 

To define the various Green's functions of a fermion in a medium, we first introduce 
Heisenberg field operators for the dynamical spinor components {a G {1,2}), 

Mr) = e^^"^•^/2$„(r)e-^^"^^/^ (14) 
^t(^) = g.p-r.+/2$t(^)e-^-^+/2. (15) 

At equal r+, these Heisenberg operators fulfill the same anticommutation relations, 
Eq. ()A11|) . as the Schrodinger operators. 

The time-ordered Green's functions of a fermion in a medium is defined in terms of the 

B 



dynamical fields only 



tGaAri,r2) = {T+Mn)^lir2)) (16) 

= A(n)V^J(r2))e(r+ - rt) - {^lir,)iPain))eirt - r+), (17) 

where the average, 

(...) = %^, (18) 

is to be taken with the appropriate ensemble, i.e. canonical or grand canonical for large 
systems or microcanonical for hadrons. This definition of the Green's function includes the 
case of zero temperature. Therefore, the conventional light-front quantization at temperature 
T = can be formulated in terms of Gq,/? as well. 

We stress that the Green's function in light-front quantization is not a Lorentz scalar, 
in contrast to equal-time quantization. In addition, for isotropic and homogeneous systems. 
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the Ga,i3 have the general form (r = ri — 



(19) 



which is completely different from the fermion propagator in equal-time field theory. While 
G is an even function of rj_, i.e. G{r~^ ,r~ ,f±) = G{r~^ ,r~ , —f±), it can be seen from the 
definition Eq. (fT^ that G is not an even function of r"*". Neither is G an even function of 
r~, because the field operators, Eqs. ()A8IA9|) . are not defined for negative k~^, so that the 
Fourier transform Ga,[3{k) has a cut along the negative axis. This cut is needed to obtain 
separate information about fermion and antifermion distributions from the Green's function. 
The analytic structure of the light-front Green's function is therefore rather complicated. 
The Green's function in momentum space is given by 



GaAk) = / d^re^'^-G^Ar) 



(20) 



and the inverse of this transformation is 



GQ,/3(r) 



(27r) 



e{k+) (e-'^^'-'/^ + e+^f'+^-fA e-'^'^-^l^+'^^-^^G^^k). (21) 



In a gauge theory, it is also necessary to include a (path ordered) gauge connector along 
the light-cone by redefining the fermion fields, 



V'a(n) ^a(ri)Pexp{ 



^9 



Pexp{ 



dr A'^{r2,r ,f2^±)} 



(22) 



(23) 



The path-ordering symbol P arranges the gauge fields A~^{r) in order of ascending r~, i.e. 
with the largest r~ to the left. This redefinition does not affect the anti-commutation 
relations Eq. (jAllll of the fermion fields. Inclusion of the gauge connector makes the Green's 
function gauge invariant in the limit r j*" — r J ^ ±0, ri _|_ — > r2,±- The time evolution of Ga,i3 
is given by 



d 



2i^G^An,r2)-{T+ Mri),P' V^^ra)) = ^..^g^^'H^ - rs) 



(24) 
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Note that a gauge invariant Hamiltonian P~ alone is not sufficient to make Ga,f3 gauge 
invariant. One also has to include the path-ordered exponentials as in Eqs. ()22|2H|1 . 

A possible subtlety arises in light-cone gauge, A'^ = 0. Even though the path-ordered 
exponentials reduce to unity, they may still influence the time evolution of Ga,i3 because 
their commutator with P~ will in general not vanish, i.e. 



=0, Pa+=o 



(25) 



In ordinary light-front quantization at T = 0, the gauge link is associated with the single-spin 
asymmetry j28|. Its role at flnite temperature is yet to be explored. 

In addition, the retarded (R) and advanced (A) Green's functions are deflned by the 
anticommutators 



^Gj^(ri,r2) 



4), 



^a(r-i),^J(r2)|)e(rj^ - 



(26) 
(27) 



These Green's functions are needed in flnite temperature perturbation theory and ohe^ 
Eq. as well. The time-ordered Green's function Ga^p is related to G^^ and G^j^ by |27 | 



GaA^) 



^ [Gl,{k) + Gi,{k)) + ^tanh (Gl,{k) - Gi,{k)) . 



(2J 



In order to obtain this formula, it is important that the Green's function can be interpreted 
in terms of transition amplitudes. This is possible, even in a gauge theory, because the 
gauge link between the fermion flelds at ri and r2 can be absorbed into a redeflnition of the 
fleld operators, see Eq. ()22l23p . Note that at T = 0, Ga^pik) = G^p{k) for positive energy 
solutions and Ga^pik) = G'^pik) for negative energy solutions. 

In the free theory, all Green's functions fulfill the integro-differential equation, 

,2" 



m 



GiT^(ri,r2)=5«,^5(^)(n-r2), 



(29) 



with boundary conditions appropriate to their physical meaning. The Green's functions 
of a fermion in an ideal gas of temperature T = 1/(3 were presented flrst in Ref. In 
momentum space, adjusted to our notation, they read 

k+ 



'^0^,13 



^2 _ ^2 _|_ ^esgn(M/c) 



(30) 
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= S^,p (P - ^sgn(«A:)vrtanh - m^)^ , (32) 

where P refers to principle value prescription. In the limit T — > 0, the light-front Feynman 
propagator for fermions turns out to be 

Spik) = 5a,p-r^ 5— — , (33) 

k'^ — m'^-^ie 

which agrees with the result obtained earlier in Ref. j22| with a different technique. 

The pole prescriptions ±zesgn(uA;) for the retarded and advanced Green's functions are 
another manifestation of the special meaning of the equal-time energy. These prescriptions 
ensure that G^°^^(r) vanishes outside the forward lightcone, while G^^^^{r) is non- vanishing 
only inside the backward lightcone. That this is indeed the case can be seen by comparing 



p (k) and ^ (k) to the well-known propagators of a scalar particle j30| , which differ only 
by the overall factor da^pk^ from the fermion Green's functions in light-front quantization. 

Most importantly, knowledge of the correct pole prescription eliminates ambiguities in 
the definition of the non-local operator 1/k'^, which appears e.g. in the free light-cone 
Hamiltonian, Eq. ()A16|) . However, the correct prescription for the l//c^-pole depends on 
the type of Green's function and on the value of the other momentum components. The 
easiest way to overcome these complications is to construct a perturbative expansion of G 
from Eq. ()24j). The 4-dimensional (5-function on the right hand side leads to 4-dimensional 
integrals in momentum space, of the kind in Eq. ()21|) . This type of perturbation theory is by 
no means identical to the usual Feynman perturbation theory in equal-time quantization, as 
can be seen immediately from the free propagators. Nevertheless, observable quantities must 

)e the same in both formulations, at least within experimental error bars. First applications 

29^ at T = indicate that this is the case. 



It is possib 
(See also Ref. 



e to introduce the chemical potential in this formalism in a covariant way. 
.) The effective Hamiltonian in thermal field theory on the light-front is 

il = u^P" - ijQ. (34) 

This operator propagates the system in proper time r. Since in equilibrium is independent 
of space time, r" = u^t. The presence of a chemical potential then modifies the generators 
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of space-time translations as if fiu" were a gauge field. The modified translation generators 
read 

= p- - ^Qu\ (35) 

They propagate the system along trajectories of constant charge and take into account that 
the medium moves as a whole as the "test particle" represented by the Green's function 
propagates from r2 to ri. The definition of the Heisenberg operators has to be changed 
accordingly, i.e. one needs to replace P~ by V~ in Eqs. p4pi5j) . 

Another remarkable property of the light-front Green's functions is, that if the theory 
is discretized on a lattice in r-space, the factor in the numerator leads to only one pair 
of fermion doublers. Naturally, this simplification comes at a price. In QCD for example, 
the light-front Hamiltonian contains non-local terms (xl/k^ and oc l/{k^Y 3l- ■'■^ order to 
calculate the Green's function on a light-front lattice, one needs to find field configurations 
with maximum statistical weight, as determined from w. This can in principle be done 
with the Metropolis algorithm j^. Unfortunately, each configuration update requires the 
evaluation of two integrals (because of the oc l/(/c+)^-terms) over the light-cone (for every 
lattice site.) It is at present not clear, if this disadvantage outweighs the absence of most 
doubler states. However, we argue that a light-front lattice is advantageous at large chem- 
ical potential, where Monte-Carlo techniques usually cannot be applied because of the sign 
problem. There is no fermion determinant in the Hamiltonian approach developed here, and 
also the chemical potential enters in a different way, see Eq. A further investigation 

of light-front lattice techniques seems therefore interesting to us. Results from transverse 
lattice calculations at finite temperature have been published recently as well js^. 

In Discretized Light-Cone Quantization DLCQ |ll8[, the alternative non-perturbative tech- 
nique in quantum field theory, it is already known that no fermion doubling problem occurs. 
Furthermore, it is possible to do calculations for a fixed value of the total charge without 
any sign problem. The output of DLCQ is the full spectrum of the theory, i.e. all masses 
of the energy eigenstates and their wavef unctions. This is of course much more information 
than is needed to calculate bulk thermodynamic properties. Nevertheless, from the DLCQ 
solution of a field theory, one can directly obtain the partition function. This was actually 
done in Ref. Q for QED(1+1). 
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Once the Green's function is known numerically from any method, one can calculate all 
thermodynamic properties of the system by integrating the relation 

In fact, it is sufficient to know G{r) for all fi and T in the limit r — 0~ to calculate the 
grand-canonical potential Q. From this, the equation of state and hence the phase structure 
of the theory immediately follow. The relation of Ga,i3 to the total charge Q will be explained 
in more detail in the next section. 

IV. THE LIGHT-FRONT DENSITY MATRIX AND GENERALIZED PARTON 

DISTRIBUTIONS 

The physical meaning of the time-ordered Green's function Ga,f3 is eludicated further by 
comparing its definition, Eq. (jl6|) . with the charge operator in Eq. ()A12|) . One observes that 
the net charge distribution in coordinate space is given by 

<5(z:) = -«G„,„(r+ ^ 0-,r|r+,r). (37) 

Summation over a is understood. The notation — > 0~ means that the limit r"*" — > is to 
be taken from below. The value of the chemical potential /i as a function of temperature T 
and total charge Q is then determined by the equation 

Q = -^ J d^rGa^^{r^ 0-,r|r+,r). (38) 

Of course, the density Q{r) is constant, if Ga,f3{ri,r2) = Ga^piji — r2). 

In this section, we introduce the light-front density matrix, which is related to the light- 
cone wavefunctions of Ref. 3| in the same way as the one-particle density matrix is related 
to the wavefunction of a system in non-relativistic quantum mechanics. If the system is a 
hadron or a nucleus, the diagonal elements of the light-front density matrix are the parton 
distributions measured in hard scattering experiments, while the off-diagonal matrix elements 



35 



This illustrates that the 



turn out to be generalized parton distributions (GPDs) j3J, 
quark and gluon distributions in a hadron are not classical densities, but quantum mechanical 
quantities. 
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Finding all light-front wavefunctions for a hadron is equivalent to solving the QCD bound 
state equation 

P -I- A/f^ 

P~\h) = -^j:^\h). (39) 

This has been accomplished for QCD in 1 + 1 dimensions 1^3], but is a formidable task in 
any 3 + 1 dimensional field theory. On the other hand, if the system is large enough, such 
as a heavy nucleus, a heavy ion collision, or a neutron star, one can hope that a statistical 
approach may be viable. In fact, for a system with a very large number of partons, say > 10^, 
the density-matrix is probably more useful a concept than the light-front wavefunctions. 

Knowledge of the density matrix enables one to calculate the expectation value of any 
single-particle operator, i. e. any operator that can be written as 

FaAr+) = J d\^i{r)Uj,{r) (40) 

in second quantization. The density matrix provides only an incomplete description of the 
system, but is sufficient for many applications of statistical physics. 

A relation similar to Eq. (j?fj) exists in condensed matter physics [2^, but in light-front 
field theory the situation is more complicated due to the presence of anti-particles. For 
example, one is tempted to interpret 

QaA'Li.L2) = -iGaA^t 0~,ri|r+,r2) = (^)^(r2)^a(ri)) |^+^^+ (41) 

as the density matrix of the net charge distribution in the system, in analogy to non- 
relativistic Fermi systems. However, since the total charge Q can have either sign, one cannot 
normalize Qa,/3{LiyL2) simply by dividing it by Q. Therefore, knowledge of Qa,i3{Li:L2) does 
not immediately enable one to calculate the expectation value of single-particle operators. 
One first needs to find separate density matrices for fermions and antifermions, qa,i3{'LiyL2) 
and g„_^(ri,r2). Then, 

^^+^^ ^ Id^nd^r2fl]l [qaALi,r2) + qaA^^i,r^)] S^^^r^-r^) ^^^^ 
with fj^l acting on r^^. 
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In the 2-component theory, the separation of fermion and antifermion distributions re- 
quires the evaluation of a Fourier-integral of the Green's function. Since Ga^i3{ri,r2) often 
depends only on the difference r = ri — r2, it is natural to introduce the variables 

R = (43) 

r = ri — r2- (44) 

The first step in finding an expression for the quark density matrix is to write down the 
operator structure of Qo,/3, 

Q^du^r^) = (45) 



|6t(A;2,A2)%,Ai)xJ,(A2)Xa(Ai 



^+ik2-(R- = )-ik^-{R+l 



\ije 



+ b^ik^, \2)d\k„ Ai)x^(A2)Xa(-Ai)e+^^^ (^-i)+^^^ (^+i) 

+ fe,A2)K^i,Ai)xJi(-A2)x«(Ai)e-^^-(^-i)-^^i-(^+5)}y (46) 

We resort to = gauge here. In order to isolate the quark distribution, one needs to 
Fourier transform over r~ . Since the k^' cannot be negative, the dS-teim disappears and 
the quark density matrix is given by 

qa,(3ik~^,R,r±) 

= ^1 dT-e^^'^'-"Q^Ari.r,) (47) 

\^ f ^^^1 e(/^+) I ^^^^ e(/;;+) 

b\k„ X,)b{k„ Ai)xJ(A2)x«(Ai)e+^fe-^^)-^+K^^.x+^^.x)-^-x5(fc+ - ^tt^) 

+ ?(fc2, A2)rf^(fci, A0xi(A2)x.(-Ai)e+^(^^+^^)-^+K^^.x-^^.x)-^^x5(fc+ - 
+ ^(^2, A2)%, Ai)x^(-A2)Xa(Ai)e-^(^^+^^)-^-K^^^.--^^.^)--5(fc+ _ ^t-]±)\\ .(48) 
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A similar expression can be obtained for the antiquark density matrix, 

qaAk^,R,f±) = ^J dr-e^^'^^-IhG^^rt - 0+, r^, r+, r^). (49) 

To obtain the correct order of the creation and annihilation operators for antifermions, the 
r"*" ^ limit of the Green's function is now taken from the other side. For definiteness, we 
shall discuss only the quark density matrix in the following. 

The density matrices qa^p (and are related to the so-called Wigner function by a 
Fourier transform over r±. The Wigner function is the quantum mechanical analog of the 
classical phase space distribution, to which it reduces in the limit h ^ In quantum 
mechanics, the Wigner function is real but not always positive. It is nevertheless a physical 
quantity that can be measured in experiment. We remark that all properties of the quantum 
mechanical density matrix, such as hermiticity and positivity of the diagonal matrix elements, 
also apply to the light-front density matrix in = gauge, since this gauge has only states 
with positive norm and no unphysical degrees of freedom. However, the Wigner function 
and the light-front density matrix have matrix elements that are off-diagonal in Fock-space. 
The object defined in Eq. (j47|) is similar to the Wigner function introduced in Ref. jsSj. 

The fermion density matrix contains all information about single quark properties. It 
depends on 6 variables and is a 2 x 2 matrix in spinor space, which can be written as a linear 
combination of Pauli spin matrices. The coefficients are the density matrices for unpolarized 
(g), longitudinal (A^g), and transverse (A^g^'^) spin distributions. 



qa,p 



ik+,R,f±) = q{k+,R,f±)6a,f3 + ALq{k+,R,f±)a3 

+ ATq\k+,R,f^)au+ATq\k+,R,f^)a2, (50) 



As already mentioned above, the diagonal matrix elements in coordinate space of 
qa,i3{k~^,R,rj_), i.e. the ones with fj_ = 0_l, are closely related to the usual PDFs 
For instance, the unpolarized collinear quark density is given by 



q(k+) = j d^Rq{k+,R,r^ = 6^). (51) 

Strictly speaking, this integral is divergent and needs to be renormalized, making q{x) scale 
dependent. At T = 0, the renormalization group equations are the QCD evolution equations. 
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FIG. 2: The off-diagonal matrix elements of the light-front density matrix are related to GPDs, 
which can he accessed in DVCS. The time-ordering shown on the left probes density matrix elements 
that are off-diagonal in momentum space, hut diagonal in Fock-space. The other time-ordering 
(right) samples matrix elements that are off-diagonal in Fock-space. 

A generalization of QCD evolution to finite temperature seems possible. In fact, in equal-time 
quantization, such a program has been started in Ref. 4^. Completely analogous expressions 
exist for Ai^rg(A;+). The parton density is normalized such that dk^q[k^) = q, the total 
number of quarks in the system. Note that 



X 




:b^{k, X)b{k,X) 



(52) 



2(27r)3 

depends only on the light-front momentum fraction x = k'^/P^, which goes to in the 
thermodynamic limit. It is therefore appropriate to consider all PDFs as functions of k^. 
The off-diagonal matrix elements of qa,i3 and are related to generalized parton distri- 



butions (GPDs) |34, 
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which can be accessed e.g. in deeply virtual Compton scatterin; 



(DVCS), see Fig. |21 The generahzed form factors H and E are defined by (see e.g. Ref. 
for details) 

dr 



An 



3^^-'^-/2(P',A;|^t(o)^(r)|P,A,) 



2P 



-u{p', a;) 



H{XX,th-' + E{X,C,t) 



+=0,rx=Oj 
I 



-a 



-A. 



2m7v 

A„Ai + ^A2 



(53) 

f/(P,Ap) (54) 



1-^ "'--^ V^^^-^) """" ^""^ ^" ^^^^ 

where P = (P'^ + P"')/2. The meaning of the four-momenta P and P' is illustrated in 
Fig. 121 The momentum transfer to the target is A*^ = P'^ — P"^ and t = A^. We use the 



-, (55) 
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kinematic variables of Ref. 



36[, i.e. X = kf /P^ and ( = 1 — P'^/P'^, and employ a frame 



with g+ = 0, where q is the four- momentum of the virtual photon. 

Since the electromagnetic current touches only one quark (at leading twist), the off- 
forward matrix element {P' , Xplip"^ {0)ip{r)\P, Xp) is diagonal in the momenta {pf,Pi.±) and 
the spin variables (Aj) of the spectator quarks. Note that the index _L refers to the coordinates 
transverse to the light-front z-direction, not to the direction of the proton. That way, the 
sum over all Pi^±, including the struck quark, yields the total transverse momentum of the 
proton, which in the final state is different from Oj^. Likewise, Aj is the projection of the 
Pauli-Lubansky spin vector onto the light-front ^-direction and therefore different from the 
quark's helicity. We remark that the usual light-front variables Xi and fcj _|_ of the spectators 
(with ^ /cj _|_ = 0_|_) are changed in DVCS, but this is a purely kinematical effect. Therefore, 
in terms of the variables Xi and ki^±, {P',Xp\ilj'^{0)ilj{r)\P,Xp) can be written as overlap of 
light-cone wavefunctions 0, |4^. In terms of the actual light-front momenta, which are 
defined here with respect to the z-direction, the DVCS amplitude is simply the light-front 
density matrix, i.e. the sum (g^,/? + ^0,13)^ momentum space. 

Loosely speaking, the density matrix in coordinate space is related to the off-forward 
matrix element (as a function of and fcg) by Fourier transformation. The precise relation 
can be read off from the phase factors in Eq. (|47p . Because of the positivity constraint on the 
light-front +-momenta, one has to distinguish four different kinematical domains j^]. By 
inspecting the (5-functions in Eq. fl47|l . one finds that for < C < X < 1, only the 6^6-terms 
contributes to the amplitude (Fig. |21 (left)), while for < X < < 1, the electromagnetic 
current pulls a quark-antiquark pair out of the proton (bd), see Fig. |21 (right). In the latter 



ar considerations 
A complete 



region, there is also a contribution of g^,^ to the DVCS amplitude. Simi 
can be made for XX < 0. This is explained in more detail in Ref. 
discussion of the relation between qa,i3,Qai3 ^"^^ GPDs is beyond the scope of this paper and 
will be published elsewhere. We shall explain only one special case here, which is particularly 
interesting. 

For C = 0, one has k~^ = kf = k2 , and the GPDs can be identified as impact parameter 
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dependent parton densities 



41( 1 ■ In the case \p = \' we find 



q{k+,b) = J dR'^6f,,^q^,^ik+,R-Xr± = 6±) (56) 



A 

q{k-',b)dk+ = J^e^'^-^-H{X,C = 0,t). (58) 

The creation and destruction operators for a quark at impact parameter b are 

6(fc+,M) = |^e-^V%+,4^,A), (59) 



6t(A;+,M) = I ±th±e+''---%\k+,h,^,X). (60) 



(2vr)2 

The relation between GPDs and fe-dependent parton densities for ( = was obtained 
previously by Burkardt j4l| from a different point of view. We stress that for C 7^ 0, GPDs 
are in general not probability distributions but density matrices, which do not need to be 
positive. Note also that the density matrix only carries information about single particle 
properties and contains no information about correlations between different partons. Some 
applications of the Wigner function in the context of DVCS have already been discussed 
in Ref. 



3a |. The light-front density matrix is a natural extension of the parton model to 
quantum mechanics: classical parton densities are replaced by a density matrix. It would be 
interesting to identify other hard processes besides DVCS which are sensitive to the quan- 
tum mechanical nature of parton distributions. Most important however is the connection 
between the density matrix and the fermion Green's function, because that establishes a 
common language for high energy scattering and statistical QCD. 



V. SUMMARY AND OUTLOOK 



We have investigated the prospects and challenges of light-front quantization in statistical 
physics and thermodynamics. Though most of the results of this paper apply to any field 
theory, we are especially interested in QCD. We have constructed the most general form of the 
statistical operator w, Eq. using only the light-front Liouville theorem and the Poincare 
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algebra. Our results generalize earlier findings jj, |7[ to also include rotations and finite 
density. In particular, we treat the chemical potential in a covariant way. Remarkably, id is 
not the exponential of the light-front Hamiltonian P~, but of the operator Ti. = UyP'^ — (iQ^ 
which propagates the system in eigentime r = Ut^r^ , i.e. along its world line. This is a 
direct consequence of cluster decomposition, i.e. of the requirement that the entropies of 
two systems separated by a large distance are additive. Since the spectrum of Ti (for (5 = 0) 
is that of the equal-time Hamiltonian in the rest frame of the system, no problems related 
to light-cone zero-modes arise in the description of phase-transitions. This finding is also 
important for the conventional light-front quantization at zero temperature: at T = 0, the 
system is in the state with the lowest eigenvalue of 7i and not in the state with lowest P~ . 

Despite all this, thermal field theory in light-front quantization is not identical to the usual 
equal-time formulation. In light-front quantization, commutation and anticommutation rules 
are defined on a light-like hypersurface and only one of the four translation generators, 
namely P~, depends on the interaction. The other three terms in u^P'^ are kinematic. It 
is therefore still possible to make use of the light-front Fock-state representation at finite 
temperature and density. The partition function of a theory for any T and especially any /i 
can therefore be obtained from DLCQ (see also Ref. j^). The output of DLCQ is the full 
spectrum of the theory, eigenvalues and eigenf unctions, which is much more information than 
the single particle properties needed to calculate thermodynamic quantities. It is therefore 
our hope to develop a reduced version of DLCQ that yields less detailed information but 
is numerically inexpensive enough to be applied in 3 -|- 1 dimensions. The basis for such a 
technique would be the light-front density matrix, which represents the incomplete quantum 
mechanical description of a system and replaces the light-front wavefunctions in statistical 
physics. The light-front density matrix is closely related to GPDs [38]. 

Even though formidable numerical challenges still exist in 3-|- 1 dimensional field theories, 
we believe that this result is important, because no other first principle technique exists to 
calculate thermodynamic properties at low T and large /i. (See however Ref. |l9| for recent 
progress in lattice QCD with finite chemical potential.) Unlike in lattice gauge theory, there 
are no problems arising from dynamical fermions in DLCQ. 

The density matrix can be constructed in an elementary way by integrating products of 
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light-front wavefunctions over the variables of all but one particle. The light-front wavefunc- 
tions can be obtained numerically by diagonalizing the DLCQ Hamiltonian. This procedure 
is formally equivalent to taking the limit r+ 0^ of the Green's function of a particle in the 
medium. The latter approach yields the most obvious generalization of GPDs to finite tem- 
perature and density. It therefore provides a common language for the description of parton 
distributions in high energy scattering and the thermal properties of a fireball. Although it 
may be difficult to measure finite temperature PDFs in deep inelastic scattering type exper- 
iments, these quantities have a profound impact on jets passing through the mediurn |lfil| 
and will also modify the fragmentation functions of heavy quarks and other particles [l7|. 

We regard the Green's functions as the fundamental quantities the theory is built upon, 
both at T = or T 7^ 0. The approach presented in this paper may therefore be regarded 
as a novel way of doing light-front quantization. The advantage of the Green's function 
formulation is that there are no ambiguities in defining the boundary conditions for the 
operator l/d~^. The pole prescription for the Green's functions are fixed by the requirement 
that positive energy solutions have to propagate into the forward light-cone, while negative 
energy solutions are to propagate into the backward light-cone. 

The fermion Green's function in the front-form has properties very different from the 
equal-time Green's function. The light-front Green's function is a 2 x 2 matrix in spinor space 
and does not transform as a Lorentz scalar. Its close relation to GPDs makes light-front field 
theory appear very similar to condensed matter physics S^]. There are, however, important 
differences between relativistic light-front field theory and nonrelativistic condensed matter 
physics: The necessity of antiparticles in a relativistic theory leads to a rather complicated 
analyticity structure of light-front Green's functions. In particular, the Green's functions 
are not defined along the negative p+-axis. This fact is important to separate particles from 
antiparticles in a theory with 2-component spinors. 

In summary, we have presented a new formalism for analyzing relativistic statistical sys- 
tems based on light-front quantization. The new formalism provides a boost-invariant gen- 
eralization of thermodynamics, and thus it has direct applicability to the QCD analysis of 
heavy ion collisions and other systems of relativistic particles. 
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APPENDIX A: LIGHT-FRONT QUANTIZATION OF THE DIRAC FIELD 

The Dirac equations are derived from the usual fermion Lagrangian L = — m)il) 
(with ip = ip'^'j^, not ?/'^7~). It is a special property of front form dynamics that the Dirac 
equations are a set of two coupled spinor equations P|, only one of which contains a time 
derivative, d~ = d/dx^ = 2d/dx~^ . This property of the light-front Dirac equation becomes 
most transparent in the following representation of the 7-matrices In 2 x 2 block matrix 
notation, 

^ /o o\ [ -2i\ [ -ta± \ 

y2t J yo J \ tax J 

where a± = (cri,(T2) are two of the three Pauli spin matrices. The operators projecting out 
the free (A+) and the constrained (A_) components of a 4-spinor take the simple form 

1 ^ /^i o\ 1 ^ [ o\ 

A+ = -7-7-^ = ' A- = 77^7- = , (A2) 

and the Dirac equations can be written in terms of 2-component spinors and 77, 

id~ip{r) = (^—ia± ■ d± — irnj rj{r) (A3) 
id'^7]{r) = (^—i(r± ■ d± + im^ i}j{r). (A4) 

The usual 4-component spinor if) is related to ip and r] by 

^{r) = f "^^"^ I . (A5) 

The free particle solutions for positive {ip{r) = u{k)e~^^^) and negative energies {ip{r) = 
v{k)e^'^^^) read, 

u{kA) = Vk^\ . I , v{k,\) = Vk^\ . I, (A6) 

k+ aA / \ fc+ A-A 
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where x\ is an eigenstates of a-^ with eigenvalue A. These spinors are normahzed 
to u{k, X)u{k, X') = —v{k, X)v{k, X') = 2m6x^x'- Note also that u\k,X)u{k,X') = 
v'^lk, X)v{k, X') = 2k^6\^x'. We also remark that Eq. ()A1|) is a chiral representation of the 
Dirac algebra, if o"3 is chosen diagonal, 

In light-front field theory, the initial conditions and the commutation relations for field 
operators are defined on a hypersurface in Minkowski space which is perpendicular to a light- 
like vector n'^. The standard choice is = (1, 0, 0, 1) (here k G {0 . . . 3}), but it is possible 
;o formulate light-front quantization in a manifestly covariant way for general light-like n'^ 
4^ . In the following, we shall use n'^ = (1,0,0, 1). The 4-velocity u'^ then still needs to 
satisfy Weldon's condition u'^n^ > ^. 

The fermion field operators for the dynamical spinor components in the Schrodinger 
picture are expanded as 

d^k 



, , , , , , (AJ 

(27r)32VA;+ ' 



A 



with r = (r~, rj_), = (A;+, /c^) and k-r = k^r~ /2 — k±- f± (see appendix iBl for a summary 
of the notations used). The creation and annihilation operators obey the anticommutation 
relations 

|6(fc,A),?(^',A')} = {d{k,X),d){k!,X')^ = {2Txf2k+5^^\k-U)h,w^ (AlO) 

so that the anticommutator of the dynamical spinor components reads 

|$„(r),$l3(r')} = 5^,^5^^\r-r:). (All) 

Remarkably, fermions are completely described on the light-front by 2-component spinors. 
This gives the theory a certain non-relativistic appearance, even though all particle and 
anti-particle states are included. 
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The formal similarity between light-front field theory and non-relativistic many-body 
physics may occasionally be misleading. For example, the operator of the conserved charge 
is 

Q = \j rfV^(r)7+^(r) = j dV$t(r)$(r) (A12) 
= E / (2v032fc+ ®^^''^ P^^-' " ' ^^^^^ 

A 

which counts the number of fermions minus anti-fermions. Here \E'(r) is the operator of the 
usual 4-component spinor field The operator '^'^'^ is not positive definite as is the case in 
non-relativistic systems. From the anticommutation relations 

|Q,$(r)} = $(r)(2Q- 1) = (2g + l)$(r) (A14) 
|g, $t(^)| = ^>\r){2Q + 1) = {2Q - l)$t(r), (A15) 

one concludes that ^!\r) creates a fermion at position r in coordinate space, while ^{r) 
destroys one at the same point. 

Furthermore, we display the expression for the light-front Hamiltonian of free fermions, 

P = d'^r^^r) ^(r) (A16) 

= E / j^^^^^^^'^'^w^ P^^^' ^) + • 

This operator is non-local because of the {id^)^^. Boundary conditions have to be chosen 
such that positive energy particles propagate only into the forward light-cone, while negative 
energy solutions propagate only into the backward light-cone. 

APPENDIX B: NOTATION 



We follow the conventions of Ref. [23], i-G. for any vector v = {v~^,v~,v±), we define 
■y+ = -y" + and v~ = — v^. The metric tensor is then off-diagonal P|, so that the scalar 
product is 

V1V2 = ^77^ + - Vi± ■ V2± = Vi+V^ + Vi_V2 - Vi± ■ V2±. (Bl) 
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Partial derivatives refer to coordinate space, 

d- 



_d_ 
dr- 
_d_ 



d 
d 



The 4-dimensional volume element is 



d^r = -dr^dr d^r\, 
2 



(B2) 
(B3) 

(B4) 



and we denote the 4-dimensional 5-function by 

5'^^\r) = 26{r+)6{r~)6{f±) = 6{r_)6{r^)6{f±) = 6{r+)6{r+)6{fi_), (B5) 

so that Jd^r6^^\r) = 1. 

In addition, it is useful to introduce a shorthand notation for spacelike vector-components, 
where spacelike means the three components r = {r~,f±) of coordinate space vectors, but 
k = {k~^, k±) for momenta. We then define 



r ■ k 



-k+ 



r± ■ ki 



d^r = dr d^r±, 

d^k = dk~^d'^k±, 

5(3) (r) = 5(r")5(rl) 

= 6{k-')6{k±) 



1 
2 
1 
2 



(27r)3 
d^r 



The Lorentz invariant phase space element is given by 



d^k 



d^k 



(27r)4 ' ' (27r)32A;- 

Furthermore, we denote the step-function by 

1 for c> 



-Q{k^ 



e(c) 



and the sign-function by 



sgn(cj 



for c < 0. 

+ 1 for c> 
for c = 
-1 for c < 0. 



(B6) 
(B7) 
(B8) 
(B9) 

(BIO) 
(Bll) 

(B12) 



(B13) 
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